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High-resolution electron energy loss spectroscopy measurements have been carried out on an op- 
timally doped cuprate Bi2Sr2CaCu208+i- The momentum-dependent linewidth and the dispersion 
of an A 1 optical phonon are obtained. Based on these data as well as the detailed knowledge of the 
electronic structure from angle-resolved photoemission spectroscopy, we develop a scheme to deter- 
mine the full structure of electron-phonon coupling for a specific phonon mode, thus providing a 
general method for directly resolving the EPC matrix element in systems with anisotropic electronic 
structures. 
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The interaction between electrons and various collec- 
tive excitations (bosons) is the central ingredient for un- 
derstanding many of the novel physical properties in 
condensed matter systems [lH3|. For simple systems 
with isotropic electronic structures, such an interaction 
(electron-boson coupling, EBC) can be characterized by 
the Eliashberg spectral function a 2 F(uj) [4j], which de- 
scribes the energy of the bosonic modes involved in 
the coupling as well as their coupling strengths. Di- 
rect experimental determination of the Eliashberg func- 
tion has proven to be important. For instance, in Pb, 
it played a crucial role in establishing the BCS the- 
ory of superconductivity j5f. Experimentally, there ex- 
ists several techniques for determining the Eliashberg 
function, such as the McMillan-Rowell inversion method 
on the tunneling data of the conventional superconduc- 
tors [5J, la], and more recently, the maximum entropy 
method (MEM) in analyzing the quasi-particle disper- 
sion kink observed in high-resolution angle-resolved pho- 
toemission spectroscopy (ARPES) measurements [7H9J. 
These approaches attempt to elucidate the EBC by prob- 
ing its effects on the electrons. On the other hand, the 
effects of the EBC on bosonic modes is not commonly 
investigated. 

For more complex systems, neither the isotropic 
Eliashberg function a 2 F(uj) nor the anisotropic Eliash- 
berg function a 2 F(oj, k) determined from ARPES is suf- 
ficient for fully characterizing EBC. A notable example 
is the high-T c cuprate, who has highly anisotropic elec- 
tronic structure and the unconventional d x i_ y 2 super- 
conductivity. By probing the electrons alone, one can 
at the best determine the a 2 F{uj, k) with the depen- 
dence on the direction k alone using ARPES and the 
aforementioned MEM analysis [lfj. However, unlike the 
case of the conventional s-wave superconductors, the so- 
determined Eliashberg function cannot be directly re- 



lated to the strength of the unconventional d x 2_ y 2 pair- 
ing. This is because the Eliashberg function determined 
from the quasi-particle dispersion and the one dictating 
the d x 2 „2 pairing strength belong to different symme- 
tries 11] . For these systems, it is necessary to resolve 
the full structure of the EBC, i.e., the matrix element 
g(fc,fc'), which characterizes the probability amplitude 
for the electron transition from k to k' induced by the 
interaction with bosons [12|, |l3|. Because the electronic 
measurements such as ARPES only contain the infor- 
mation of the electron final state fe', with the contribu- 
tions from the different initial states k integrated, it is 
in general impossible to determine the full structure of 
g(k, k') by probing electrons alone. The manifestation of 
the EBC on the bosonic modes needs to be probed. 

In this paper, we develop a method to resolve the 
full structure of EBC by probing the bosonic modes 
with high resolution electron energy loss spectroscopy 
(HREELS), in combination with the existing ARPES 
data for the electronic properties. The momentum- 
dependent linewidth as well as the dispersion of an A\ 
optical phonon was recorded as a function of the direction 
relative to the crystalline orientation. The initial and fi- 
nal states involved in the coupling and the corresponding 
strength are simultaneously determined. The developed 
scheme is general and applicable to many systems with 
anisotropic electronic structures. 

We have used a prototype high-temperature super- 
conductor (HTSC) Bi2Sr2CaCu208+«5 (referred to as 
BSCCO in the following) 0, [lj| as a test system. The 
electron band structure, Fermi surface, and surface prop- 
erties have been well characterized [16| . Many experi- 
mental studies have revealed the signature of the inter- 
action between electrons and a collective bosonic mode. 
ARPES shows the electron self-energy renormalization in 
the form of kinks in the dispersion at the energy scales of 



~70 meV below the Fermi energy in the nodal re gion 1171- 
|20( and ~40 meV near the antinodal region 13, l2l| . 
Direct measurements of phonon spectra have been car- 
ried out with Raman spectroscopy |22J ]. infrared spec- 
troscopy (IR 1 1231 . HREELS [24j, inelastic neutron scat- 
tering (INS) |25[ and inelastic x-ray scattering (IXS) [26j 
for BSCCO and other cuprates HTSCs. However, a 
systematic study of the momentum-dependent phonon 
renormalization is still missing and it is needed to eluci- 
date the details of EBC. 

The experiments were carried out in a commercial sys- 
tem containing variable temperature scanning tunneling 
microscope (STM) and HREELS (LK-5000). The base 
pressure was better than lxlO~ 10 mbar. The supercon- 
ducting property of the optimal doped BSCCO single 
crystal was carefully characterized with the transition 
temperature at T c = 91 K. The sample was cleaved 
in situ in vacuum and the surface was characterized 
with STM and low energy electron diffraction (LEED). 
Fig. 1 (a) shows that the cleaved surface was atomically 
flat with the Bi-0 plane exposed. The cleavage quality 
was further verified by the presence of the (1x5) super- 
structure that was visible in both the LEED patterns 
and STM images, which also indicates the relative angle 
of the crystalline orientation to the measurements geom- 
etry. The HREELS measurements were taken with the 
electron incident angle of 65°, and excitation energy var- 
ied from 3.5 eV to 50 eV to achieve different detecting 
depth. All the observed spectral features showed no exci- 
tation energy dependence for a fixed momentum transfer 
parallel to the surface. Here we report the results taken 
with the excitation energy of 50 eV offering a wide mo- 
mentum transfer range. The HREELS sample stage was 
coupled to a constant-flow liquid helium cryostat. And 
the measurements were taken at different temperatures 
from 60 to 300 K, across the superconducting T c . 

Angle-resolved HREELS of the optimal doped BSCCO 
sample were measured along different directions relative 
to the crystalline orientation, as shown in Fig. 1. At Bril- 
louin zone (BZ) center (the in-plane momentum transfer 
q = 0), two main features are resolved near 50 meV (Fi) 
and 80 meV (F 2 ), respectively. Feature Fi corresponds 
to the out-of-plane vibration of oxygen atoms in Bi-0 
plane 22j-|24j, with two shoulders at both sides with en- 
-41 meV and ^62 meV that have been reported 
Feature F 2 , to be focused on in the follow- 
ing, appears as a main peak at ^80 meV with a broad 
tail centered at ^88 meV. The vibration along c axis of 
the apical oxygen atoms in the Cu-0 semi-octahedra is 
responsible for the main peak 22j-|24j, while the high- 
energy tail might be related to the disorder [22]. Also 
considering the selection rule of HREELS, we are able to 
rule out other energy loss mechanisms that fall into this 
energy range, such as the phonons of in-plane breath- 
ing modes of oxygen in CU-O2 plane since they are not 
dipole active along the surface normal and therefore un- 
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FIG. 1: (color online) (a) STM image of the cleaved BSCCO 
surface. The left inset shows the 1x5 superstructure with 
the " x 5" direction aligning to [110] orientation. The right 
inset shows the LEED patterns taken with beam energy of 
61 eV at room temperature. The line profile along AB is 
displayed in the lower panel, which indicates a step of 1.5 nm 
corresponding to the interspacing between two cleavable BiO 
planes, (b) Schematic drawing of the scattering geometry for 
the angle-resolved HREELS measurements, (c) and (d) The 
HREELS spectra measured with ip — 45° and 0°, respectively, 
with the excitation energy of 50 eV at 60 K. The bars (red) 
are guides to the eye for the phonon energy shift. 



detectable by HREELS at BZ center. Along any mea- 
surement direction, we detect no obvious temperature 
dependence across superconducting T c . 

As shown in Fig. 1 (c), the most prominent character- 
istic of the momentum-resolved HREELS presents along 
nodal direction (ip = 45°) where the energy dispersion 
of the 80-meV feature F 2 softens significantly from BZ 
center towards the boundary. Such a behavior has been 
quantitatively analyzed 27] . The momentum-dependent 
phonon energy and linewidth (full-width at half maxi- 
mum, FWHM) along different directions relative to the 
crystalline orientation are shown in Fig. 2. An abrupt 
decrease from 81 meV to 74 meV occurs to the phonon 
energy with an onset at qn ~ 0.45 A -1 (0.54 7r/a with a = 
3.82 A). With the same onset, the FWHM increases sud- 
denly from 8 meV at q|| ~ 0.45 A -1 to 18 meV at q\\ ~ 0.8 
A~\ indicating the intrinsic strong interaction between 
the apical oxygen phonon and electrons in BSCCO. 

The phonon dispersion and linewidth data show dra- 
matic dependence on the measurement direction with re- 
spect to crystalline orientation (ip). As ip decreases from 
45° (the nodal region), the anomalous behavior persists 
until ip gets close to the antinodal region (tp < 10°), when 
neither the energy nor the linewidth shows prominent 
momentum-dependence at temperatures above or below 
T c [also see Fig. 1 (d)]. Therefore the phonon disper- 
sion and linewidth measured with different ip values can 
be divided into two groups, respectively. Near the nodal 
region, the phonon energy decreases and the linewidth 
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FIG. 2: (color online) (a) The energy dispersion of the apical 
oxygen phonon at 60 K. The different azimuth angels of qy 
are illustrated in the upper inset in k space relative to the 
Fermi surface. When ip >10°, there is a sudden drop with 
the onset of qy ~ 0.45 A -1 . The lower inset compares data 
taken at different temperatures, indicating the same energy 
softening across superconducting T c . (b) The dependence of 
linewidth on qy along different azimuth angels. There is a 
sudden broadening when ip >10°, with the onset aligning to 
that of energy dispersion at qy ~ 0.45 A -1 , as indicated with 
the black arrows. Similarly, the lower inset shows the un- 
changed linewidth dependence on qy across superconducting 
T c . All the dotted lines are guides to the eye. Due to the ex- 
istence of the (1x5) superstructure, there are two irreducible 
nodal directions (p — ±45°). Our measurements show iden- 
tical momentum-dependence along both directions. 



increases suddenly with the same onset. In contrast, 
near the antinodal region, the phonon linewidth increases 
slowly by less than 3 mcV while the energy are almost 
momentum-independent. 

Such a dramatic behavior seen in the phonon spectra 
must reflect the underlying structures of the electron- 
phonon coupling (EPC) for this mode. In general, 
the momentum-dependent broadening of the phonon 
linewidth 7epc(q) induced by EPC can be connected to 
the electron polarizability x(q, &q) by the EPC matrix 
element g(q) [J]: 



7epc(<?) = -\g(q)\ 2 Im[x{q,n q )}, 



(1) 



where Q q is the phonon energy. And we have assumed 
g(k, k') = g(k' — fc) = g(q) 1 which is a good approxima- 
tion for the particular out-of-plane apical oxygen phonon 
mode. The imaginary part of x(q, £l q ) is given by: 



Im[ X (q,n q )]= [ 
Jbz 



dfc 



dw 

2^ 



A(k, w)A(k + q,co + Cl q ) 



„ f dfc y[f(e k+q )-f(e k )) 



where A(k, u>) is the spectral function of quasi-electron, 
f(ui) is the Fermi distribution function, and e k is the en- 
ergy dispersion of the quasi-electrons. We have further 
simplified the expression in the second line of Eq. 2 by 
assuming that the spectral function of the quasi-electron 
has the form of A(k,u) « (r)/2)/[(w - e fc ) 2 + (?7/2) 2 ], 
with a parameter r\ characterizing the life-time of the 
quasi-electron due to the many body interactions. The 
quasi-particle dispersion efc of the BSCCO system has 
been measured previously (ARPES) and fitted in a tight- 
binding phenomenological model [lg, |28[. Im[x(q, Cl q )] 
can then be calculated numerically. Fig. 3 (a) shows 
the calculated Im[x(g, fi q )], which displays rather com- 
plicated features originated from the highly anisotropic 
electronic structure. 

In addition to jepc, the experimentally observed 
linewidth also include a momentum-independent back- 
ground 70, which is most likely originated from surface 
roughness of the cleaved sample and the instrumentation 
broadening: 



7exp(<?) =70 +7EPC(<?)- 



(3) 



We estimate such a background using the value of the 
experimental data near q = 0: 7o = 7cxp(<Z — > 0), 
since there should have been no EPC without momen- 
tum transfer. With the background subtracted, the EPC 
matrix element can then be determined straightforwardly 
with Eq. 1-3. We note that only the electron transition 
satisfying both energy and momentum conservations, i.e., 
fc' — fc = q and e k ' — £fe = ti&, q , contributes to the EPC- 
induced broadening. This imposes a stringent constraint 
on the possible initial and final states of electron for 
a given phonon momentum q, as shown in Fig. 3 (c). 
Therefore, the determined |g(q)| 2 should be considered 
as a subset of the more general \g(k, k')\ 2 . 

The determined EPC matrix element |.g(q)| 2 are shown 
in Fig. 4 (a) along different directions. With the least 
square method, the data can be well fitted by: 



\g q \ 2 = [(2a+b)-a(cosq x +cosq y )-bcosq x cosq y ] 2 , (4) 

which is consistent with an A\ phonon mode that cou- 
ples to electrons in the lattice with C^ v symmetry. Here 
a and b are constants, representing the coupling strength 
between an atom displacement and its nearest and next- 
nearest neighboring sites, respectively. It can be seen 
that the experimental data, in particular, the dramati- 
cally different behaviors for nodal and anti-nodal direc- 
tions, are quantitatively reproduced by our model. 

To further test our analysis, we calculate the phonon 
softening induced by EPC: 



(2) Re[ X ] 



Afi(qr) = -\g q \ 2 Re[ X (q,n q )} 
dfc (efc - e k+q 



(5) 

m q ) [f(e k+q ) - /( Cfc )] 



BZ 



(2^ 



(efc — Cfc+q 



Ml q ) 2 
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FIG. 3: (color online) (a) and (b) The imaginary and real 
parts of x(q,£l q ) calculated with Eq. 2 at 60 K with n = 10 
meV (not a sensitive parameter) for different <p, respectively, 
(c) The energy contour plot (dotted lines) of the quasi-particle 
band near Ef of BSCCO. The initial (brown areas) states 
contributing to the EPC-induced broadening for qy ~ 0.47 
A -1 at tp = 45° can be determined by superimposing the 
density plot of the integrating function in Eq. 2 onto (c). 
Connecting by q (indicated by the arrows), the final states 
(pink areas) correspond to an energy gain that equals to Ml q . 
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FIG. 4: (color online) (a) \g q \ 2 determined by Eq. 1-3 and the 



fitting (solid lines) to Eq. 4 with a = 0.6668 and b = 0.3949. 
The inset shows the measured phonon linewidth as well as the 
data calculated with \g g \ 2 for different tp. (b) The measured 
phonon energy with the data calculated with Eq. 5. Here 
we have assume a bare phonon dispersion f2 barc = 83.4 — 
1.65(cos q x +cos q y ) — 0.51 cos q x cos % + 0.55 (cos 2g;E+cos 2q y ), 
as plotted in broken lines. 



where Re[x(q)] is the real part of electron polarizability, 
as plotted in Fig. 3 (b). Figure 4 (b) shows the experi- 
mentally observed phonon dispersion that can be repro- 
duced by assuming a bare phonon dispersion f2 bare (q) 
as the experimental background. In particular, the 
anisotropy of the phonon softening is reproduced: when 
tp ~ 45°, the phonon energy softens dramatically at the 
onset qii ~ 0.4 A -1 , while such a softening gets much 
weaker when tp ~ 0°. In Fig. 4 (b), Q e (g) represents 
the phonon dispersion if the EPC is turned off in our 



model, which is almost isotropic. 

In conclusion, the measured phonon structure com- 
bined with a detailed knowledge of the electronic struc- 
ture from ARPES can be used to determine the EPC ma- 
trix element \g(q)\ 2 directly. Our approach is completely 
general and could be applicable for other complex sys- 
tems with highly anisotropic electronic structures. Our 
study also highlight the necessity of probing bosons for 
revealing the full structure of EBC in these complex sys- 
tems. 
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